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Abstract. In this paper we semianalyticaly evaluate in- 
fluence of the exponential decay of magnetic field on the 
fate of isolated neutron stars. The fact of ROSAT obser- 
vations of several X-ray sources, which can be accreting 
old isolated neutron stars gives us an opportunity to put 
some limits on the parameters of the exponential decay. 

We argue, that, if most part of neutron stars have ap- 
proximately the same decay and initial parameters, then 
the combinations of the bottom magnetic momentum, fib, 
in the range 10^^ — 10^^'^ G cm^ and characteristic time 
scale, td, in the range ~ 10^ — lO^yrs for standard ini- 
tial magnetic momentum, /ip — 10'^*'Gcm'^, can be ex- 
cluded, because for that sets of parameters neutron stars 
never come to the stage when accretion of the interstellar 
medium on their surfaces is possible even for low velocity 
of neutron stars and relatively high density of the interstel- 
lar medium. The region of excluded parameters increases 
with fiQ decreasing. 
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1. Introduction 

Evolution of neutron stars (NSs) can be called "magneto- 
rotational" , because all main astrophysical manifestations 
of these objects are determined by their periods and mag- 
netic fields. 

Four main regimes, as described for example in 
Lipunov (1992), are possible for isolated NSs: ejector, 
when a star represents a radio pulsar, or a dead pul- 
sar, spinning down due to magneto-dipole radiation; pro- 
peller, when surrounding captured matter cannot pene- 
trate through the centrifugal barrier; accretor, when mat- 
ter can reach the surface, and the NSs appear as an X-ray 
source; and georotator, when gravitation becomes insignif- 
icant, because magnetic pressure dominates everywhere 
over the gravitational pull, and geo-like magnetosphere is 
formed. 



Magnetic field decay in NSs is an uncertain subject. 
Many models were suggested (see, for example. Ding et 
al, 1993; Jahan Miri & Bhattacharya, 1994 ). The only 
strong observational result is that for the exponential (or 
nearly exponential) decay characteristic time scale, td, is 
longer than ~ 10^ yrs, because no expected effects of decay 
are observed in radio pulsars (Lyne et al. 1998). 

Field decay was used in the case of old accreting iso- 
lated NSs by Konenkov & Popov (1997) and Wang (1997) 
to explain properties of the source RX J0720-3125. And it 
seems, that if this source really is an accreting old NS and 
if it was born as a normal radio pulsar (with small period, 
<C 1 s, and high magnetic field, ~ lO^'^ G), its properties 
can be explained only if decay is working (or if this source 
was born with such unusual parameters) . Recently the in- 
fiuence of the field decay in isolated NSs was investigated 
in Golpi et al. (1998) and Livio et al. (1998). An attempt 
to include field decay into population synthesis of isolated 
NSs was made in Popov et al. (1999). 

Here we try to put some limits on the parameters of the 
exponential field decay, assuming, that old isolated NSs 
are observed as accreting X-ray sources (Walter, Wolk & 
Neiihauser 1996; Haberl et al. 1996, 1998; Neiihauser & 
Triimper 1999; Schwope et al. 1999), neglecting the pos- 
sibility, that all of them can be highly magnetized NSs, 
"magnetars" (see a brief discussion on this assumption in 
Popov et al., 1999). 

2. Calculations and results 

The main idea of our work is to calculate the ejector time, 
Ie, i.e. a time interval spent by a NS on the ejector stage, 
for different parameters of the field decay and standard 
assumptions on the initial parameters of a NS, and to 
compare that time with the Hubble time, tn- 

For constant field ts monotonically depends upon NS's 
velocity and ISM density: 



tEi^J' — const) ~ 10^/i3Q^n ^^^wioyrs. 



(1) 
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li tE for decaying field for some sets of parameters is 
greater than tn ~ 10"'^° yrs even for relatively high con- 
centration of interstellar medium (ISM), lcm~^, and 
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Fig. 1. Tracks on P-B diagram. Tracks are plotted for 
bottom polar magnetic field 8 • 10^° G, initial polar field 
2 • 10^2 Q velocity lOkms'S ISM density 1 cm'^ and 
different td- The final points for different tracks are cor- 
responding to the following moments: for tj, = 10'' yrs to 
~ 10^" yrs; for td = 10^ yrs to ~ 10^° yrs; for td = 10^ yrs 
to - 1.5 • 10^ yrs; for td = 10^" yrs to - 2 ■ 10^ yrs. Here 
we used po = 0.020 s. The line with diamonds shows the 
ejector period, Pe- 



small spatial velocity of a NS, ~ lOkms^^ (this velocity 
is about the sound velocity in the ISM), than for that pa- 
rameters of the decay, td and /ib, no accreting NSs would 
be observed, and such sets can be excluded as progenitors 
of old accreting isolated NSs. In addition, if we assume, 
that most part of NSs evolve with similar td and /ih and 
are born with some typical initial period, po, and initial 
magnetic momentum, /xo (that is supported by radio pul- 
sars observations), then such sets of parameters should be 
totally excluded from models of the magnetic field decay. 

As the first approximation only ts can be taken into 
account, because the time interval, spent on the propeller 
stage, tp, is uncertain, but usually shorter, than ts (see 
Lipunov & Popov, 1995). These authors argue, that for 
non-decaying magnetic field tE is always longer than tp. 
For decaying magnetic field the situation can be different, 
and tp can be comparable to f^- So we obtain here "upper 
limits" , i.e. we calculate very strong limits on parameters 
td and /Ltf), which can become only wider if one takes into 
account the propeller stage. 

We assume exponential field decay: 

At = /io • e"*/*^ > Mf, (2) 

where /ig is the initial magnetic momentum (/i = 
^BpR^g, here Bp - polar magnetic field, and Rns~ NS 
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Bottom magnetic momentum 

Fig. 2. Ejector time, tE, (in billion years) vs. bottom 
magnetic momentum. Curves are shown for two values 
of the initial momentum: 10^°G cm^ (upper curves) and 
10=^1 Gcm^. 

radius), td - characteristic time scale of the decay, and /U(, 
- bottom magnetic momentum, which is reached in: 

tcr = id • In (^^) . (3) 

After that moment magnetic field is assumed to be con- 
stant. 

In Fig. 1 we show, as an illustration, evolutionary 
tracks of NSs on P — B-diagram for v = lOkms"^ and 
n = 1 cm~^. Tracks start at t = 0, when p = po = 0.020s 
and n = muo = lO'^^Gcm'^, and end a,t t = tH = 10^° yrs 
(for td = 10^°yrs, td = 10'^ yrs and for constant magnetic 
field) or at the moment, whenp = pe (for td = 10^ yrs and 
td = 10^ yrs). The line with diamonds shows p = pe{B). 

As far as the accretion rate from the ISM is small (even 
for our parameters), less than ^ lO'^^gs"'^, no influence 
of accretion onto decay was taken into account (see Urpin 
et al., 1996). 

The ejector stage lasts until the critical ejector period, 
Pe, is reached: 

. 1/2 _i/4 1/2 

PE = 11.5 ^3^, n ''^^ s, (4) 



where vio = + tif/lOkms ^. Vp- spatial velocity of 

a NS. Here the sound velocity, Vg, was taken into ac- 
count, but as far as normally NSs spatial velocities arc 
higher than lOkms"^ and the sound velocity outside hot 
low density ISM regions is lower than lOkms"^, we have 

^Vp + Vg ^ Vp. n is a. concentration of the ISM. 

Initial period should be taken to be much smaller than 
Pe- We used po = s. Variations of po, if it stays much 
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less than pE, have Uttle influence on our results, i.e. in 
that case Ie is determined only by pe and history of the 
decay. We calculated spin-down according to magneto- 
dipole formula (but other regimes are possible, see Beskin 
ct al. 1993 for a review): 



dp _ 2 47rV^ 
di 3 plc^ ' 



(5) 



where fj, can be a function of time. 

For our estimates we assumed constant velocity of NSs, 
V, equal to lOkms"^ and constant ISM concentration, n, 
equal to 1 cm~^. These conditions give us a lower limit on 
ts, because normally velocity is significantly higher (frac- 
tion of slow velocity NSs is less than few percents, see 
Popov et al., 1999), and ISM density is smaller than the 
specified values (partly because high velocity NSs spend 
most of their lives in low density regions far from the 
Galactic plane). 

After simple algebra one can obtain a formula for tE, 
depending upon ta, /xq, v, n and ^b' 



-td ■ In 



tE 




+ (l 



, tE < tcr 
^-ItcrlU^ , tE 



where coefiicient A is determined by the formula: 



A = 



3Ic 



2ij,oV2vM 



-1^17 r -1 -l/2,V-l/2 



s, 



(7) 



where M can be formally determined as a combination of 

intrinsic NS's parameters, its velocity and ISM concentra- 
tion using the Bondi formula even if a NS is not on the 
accretor stage: 



M 



(8) 



Results of tE calculations for several values of /io and 
td are shown in Fig. 2. Horizontal regions in the left parts 
of all curves appear because there a NS reaches the critical 
period p = pe before it reaches bottom magnetic momen- 
tum &t t = tcr, and tE is not depended upon the rest field 
decay history. On the right side from the maximum on 
the curves a NS reaches fi = fib with period significantly 
less than PE{tib), and it takes a long time after t = tcr to 
reach p = pE (this time is longer for lower Hb). On the left 
side from the maximum, but before the horizontal region, 
a NS reaches fj, = fib with a period close to pe, and it is 
closer if fib is lower, so very quickly after t = tcr & NS can 
spin down to pe ■ The first point at the right corresponds 
to the bottom momentum equal to initial momentum, i.e. 
to the case without decay. 

One can see from that figure, that for some combi- 
nation of parameters tE is longer than the Hubble time. 




10" 10" 
Bottom magnetic momentum 

Fig. 3. Characteristic time scale of the magnetic field 
decay, td, vs. bottom magnetic momentum, fib- In the 
filled region tE is greater than 10^'^yrs. Dashed line corre- 
sponds to tu = td ■ ln{fio/fib), where tn = 10^" years. 
Solid line corresponds to PEilJ-b) = pit = tcr), where 
*d-ln {fiol fib) - Both lines and filled region are plotted 
for fif) = lO'^'^Gcm^''. Dot-dashed line is the same as the 
dashed one, but for = 5 • 10^^ Gcm^. Dotted line is a 
border of the 'forbidden" region for fiQ = 5 ■ 10^^ Gcm^. 
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Fig. 4. Characteristic time scale of the magnetic field 
decay, td, vs. bottom magnetic momentum, fib- In the 



filled region tE is greater than 10 '^yrs. Dashed line cor- 



responds to tn = td ■ ln(/io/M&)j where tn = W^'^ yrs. 
Solid line corresponds to pE^fib) = p{t = tcr), where 
tcr = td • \n{fio/fib)- Both lines and region are plotted 
for fio = 1029Gcm-3. 
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It means, that NSs never evolve further than the ejector 
stage. 

We argue, that as far as accreting isolated NSs are 

observed, combinations of td and fii, for which no accreting 
isolated NS appear can be excluded. We plotted it in Figs. 

3 and 4. 

Filled regions represent space of the parameters where 
ts is longer than lO^'^yrs, so in that region a NS never 
come to the accretor stage, and doesn't appear as accret- 
ing X-ray source. With the fact of observations of accreting 
old isolated NSs by ROSAT this region can be called "for- 
bidden" for selected parameters of the exponential field 
decay (and for specified /io). 

In the "forbidden" region in Fig. 3, which is plotted 
for fiQ = 10^°Gcm^, all NSs reach the bottom field in a 
Hubble time or faster, and evolution on the late stages of 
their lives goes on with the field equal to the bottom. The 
left side is determined approximately by the condition: 

PEiPb) =p{t = tcr)- (9) 

Small difference betwcien the line, correspondent to the 
condition above, and the left side of the "forbidden" region 
appears because a NS can slightly change its period even 
with the momentum /i = /ifc, but due to small value of the 
field angular momentum losses are also small. 

The right side of the region is roughly determined by 
the value of /it,, with which a NS can reac;li the ejector stage 
with any td, i.e. this /if, corresponds to the minimum value 
of IJ.0 with which a NS reach the ejector stage without field 
decay. 

NSs to the right from the "forbidden" region leave the 
ejector stage, because their field cannot decay down to 
low values, and spin-down is fast enough during all their 
lives as ejectors, because the bottom magnetic momentum 
there is relatively high. To the left from the "forbidden" 
region the situation is different. Spin-down of NSs is very 
small and they leave the ejector stage not because of spin- 
down, but due to decreasing of p^, which depends upon 
the magnetic momentum. 

Dashed line in Fig. 3 shows, that for all interesting 
parameters a NS with /xq — 10"^° G cm"^ reach /it in less 
than 10^" yrs. Dot-dashed line shows the same for /Uq = 
0.5- 10'''°Gcm^\ 

On Fig. 3 we also show the "forbidden" region and the 
line of reaching /if, for /io = 0.5 • 10''" Gcm^. 

Fig. 4 is plotted for /io = lO^^Gcm^. For long td, > 

4 • 10^ yrs, a NS again is not able to leave the ejector 
stage. It happens because the magnetic momentum can't 
decrease down to small value of /x (nearly /if,), and pe is 
not decreasing enough. 

3. Discussion and conclusions 

We tried to evaluate the region of parameters, forbidden 
for models of the exponential magnetic field decay in NSs 



using the fact of observations of old accreting isolated NSs 
in X-rays. 

If the main fraction of NSs have nearly the same initial 

parameters and parameters of the decay, then the inter- 
mediate values of td (~ 10^ — 10* yrs) in combination with 
the intermediate values of /if, (~ 10^* - 10^^-^ G cm^) for 
/io ~ IO"'" Gcm^ can be excluded, because for that sets of 
parameters NSs spend all their lives on the ejector stage, 
never coming to the accretor stage. 

As one can sec in Fig. 2, for higher /io NSs should reach 
ts even for td < 10* yrs, for smaller - the "forbidden" re- 
gion should become wider. Results are depended on the 
initial magnetic field, /io, ISM concentration, n, and NS 
velocity, v, so one can say, that the observed accreting iso- 
lated NSs come from, for example, the objects with high 
initial magnetic field, and the rest are not visible because 
they are in the forbidden region. To explore this idea in 
details population synthesis of NSs for realistic distribu- 
tions of V, /io and n is needed. But we can say immediately, 
that the idea of obtaining accreting old isolated NSs from 
initially high field objects is not very promising, because 
the fraction of high field NSs can not be large (basing on 
radio pulsars observations), and as far as the fraction of 
low velocity NSs is not more than several percents (Popov 
et al. 1999) and the volume fraction filled with relatively 
high density ISM is also small, accreting old isolated NSs 
should come from "typical" population, i.e. from NSs with 
/io about 10^° G cm^ or less. 

Actually, limits that we obtained are even stronger 
than they are in nature, because we didn't take into ac- 
count, that some significant time (in the case of field de- 
cay) a NS can spend on the propeller stage (spin-down 
rate at this stage is very uncertain, see the list of formu- 
lae, for example, in Lipunov & Popov 1995 or Lipunov 
1992). Calculations of this effect, and calculations for dif- 
ferent models of non-exponential field decay are subjects 
for future work. 

We cannot say anything about parameters of field de- 
cay in the case of accretion in close binaries, because there 
situation is completely different, and our results cannot 
be applied to millisecond radio pulsars or other objects in 
close binary systems. 

We note, that there is another reason due to which 
very fast decay down to small values of /if, also can be 
excluded, because it leads to huge amount of accreting 
isolated NSs. This situation is similar to "turning-off" the 
magnetic field of a NS (i.e., quenching any magnetospheric 
effect on the accreting matter), and for any velocity and 
density distributions we should expect significantly more 
accreting isolated NS than we have from ROSAT observa- 
tions (of course, for high velocities X-ray sources will be 
very dim, but close NSs can be observed even for velocities 
- 100 km s-^). 

So, the existence of several old isolated accreting NSs, 
observed by ROSAT (if it is the correct interpretation of 
observations), can put important limits on the models of 
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the magnetic field decay for isolated (without influence 
of accretion, which can stimulate field decay) NSs, and 
models, from their side, should explain the fact of obser- 
vations of ^ 10 accreting isolated NSs in the solar vicinity. 
We cannot discuss numerous details of connection between 
decay parameters and X-ray observations of isolated NSs 
without detailed calculations, we just tried to show, that 
this connection should be taken into account and made 
some illustrations of it, and future investigations in that 
field are wanted. 
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